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Gravitational lensing of the Cosmic Microwave Background (CMB) encodes cosmological infor- 
mation in the observed anisotropics of temperature and polarization. Accurate extraction of this 
additional information requires precise modeling of the covariance matrix of the power spectra of 
observed CMB fields. We introduce a new analytical model to describe the non-Gaussian structure 
of this covariance matrix and display the importance of second-order terms that were previously 
neglected. When compared with direct numerical simulations our model captures parameter errors 
to better than a few percent for cases where the non-Gaussianity causes an order unity degradation 
in errors. We also provide a detailed comparison between the information content of lensed CMB 
power spectra and ideal reconstruction of the lensing potential. We illustrate the impact of the 
non-Gaussian terms in the power spectrum covariance by providing Fisher errors on the sum of the 
masses of the neutrinos, the dark energy equation of state, and the curvature of the Universe. 

PACS numbers: 98.62. Sb, 98.70.Vc, 95.36.+X 



I. INTRODUCTION 

Gravitational potentials of large-scale structure gen- 
erate a deflection of the trajectories of the Cosmic Mi- 
crowave Background (CMB) photons, an effect known as 
CMB lensing [T]-[3] (see [4] for a review). After its initial 
detection in cross-correlation with large-scale structure 
0, CMB lensing has now been detected with high 
significance in high-resolution observations from the At- 
acama Cosmology Telescope [71 [8] and the South Pole 
Telescope [9 . 

CMB lensing generates a characteristic statistical sig- 
nature that makes the CMB sensitive to cosmological pa- 
rameters which directly influence the growth of cosmic 
structure. This breaks the angular diameter degeneracy 
in the unlensed CMB and improves constraints on param- 
eters such as neutrino masses, the dark energy equation 
of state, and the curvature of the Universe p!QHT4] . 

Mathematically, CMB lensing is described as follows. 
We introduce a vector field d(n) (the deflection field) 
such that the lensed temperature T(n) and unlensed tem- 
perature T(n) are related by 

T(ii) = f(n + d(n)) (1) 

and analogously for the Stokes parameters Q(n), U{rv) 
which describe linear CMB polarization. To lowest order 
in perturbation theory, the deflection field d(n) is the 
gradient of a scalar lensing potential (i.e. d(n) = V^(n)) 
which can be written as a line-of-sight integral: 

^(n) = -2 / ^r^ ^f^77^. ^(xfi,^), (2) 

J X{Vrec)X{V) 
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where ^ is the Newtonian potential, 77 is conformal time, 
77rec is the epoch of last scattering, and x is the angular 
diameter distance in comoving coordinates. 

CMB lensing modifies the Gaussian statistics of the 
unlensed CMB by generating a correlation between the 
primary field and its gradient [15]. It also modifies the 
shape of the temperature and ^-mode polarization power 
spectra, and generates a nonzero B-mode power spec- 
trum. This leads to two different statistical techniques 
for extracting cosmological information from CMB lens- 
ing. First, we can simply make precise measurements of 
CMB power spectra (especially the B-mode power spec- 
trum), which will include lensing contributions. Second, 
we can reconstruct the lensing potential (j) using correla- 
tions between the primary and its gradient, providing a 
new cosmological observable [T6H2Q]. 

Accurate analysis of CMB anisotropics requires correct 
modeling of the covariance matrix of the lensed power 
spectra. The lensed CMB is not a Gaussian field, and so 
its power spectrum covariance is nontrivial; in particular 
the off-diagonal correlations are important. Calculations 
of the non- Gaussian covariance of the lensed power spec- 
tra have been performed in both flat sky [211 [22] and 
full sky [23^ cases, but these calculations make the ap- 
proximation that some high-order terms in the lensing 
potential are negligible. 

In the advent of low-noise and high-resolution CMB ex- 
periments that will be able to probe polarization of the 
CMB at the arcminute scale (SPTPol, ACTPol [24], PO- 
LARBEAR [25^), it becomes necessary to assess the va- 
lidity of the current approximations for the non- Gaussian 
power spectrum covariance, and study the impact on cos- 
mological parameter estimation. The purpose of this 
paper is therefore to investigate in detail the impact 
of the non-Gaussianities induced by CMB lensing and 
quantify the information contained in the power spectra. 
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We introduce a new semi-analytical approach to compute 
the power spectrum covariance matrix and validate our 
model by Monte Carlo simulations. 

In Sec. |n[ we describe the simulations we performed 
to estimate the power spectrum covariance matrix and 
present our semi-analytical model. In Sec. |III[ we in- 
troduce a model independent way of characterizing the 
relative information content of lensed power spectra and 
idealized direct reconstruction. Finally, in Sec. \P7\ we 
apply this characterization to specific cosmological model 
parameters. 

Throughout the paper, we use a fiducial fiat ACDM 
model with the following parameters: 

{n^h^ n^h^, h, r, n„ 10^ A„ m,} (3) 
= {0.1096, 0.0226, 0.693, 0.089, 0.964, 2.419, 0.58 eV}. 

We chose a large fiducial value of rrijy so as to be testable 
with CMB lensing in the near future. 



II. NON-GAUSSIAN POWER SPECTRUM 
COVARIANCE 

To characterize the information content in CMB power 
spectra we first require an accurate characterization of 
the covariance matrix between the band power estimates 
of the temperature and polarization fields. Previous ana- 
lytic characterizations have been based on a perturbative 
expansion of the effect of lensing [211 ES] which breaks 
down in the damping tail. In this section, we first con- 
duct a suite of simulations to characterize the covariance 
and then develop analytic tools which characterize its 
main features. 



A. Simulations 

We simulate the lensed CMB on the full sky using 
the following procedure. We first make a realization 
of the unlensed CMB fields T, E and lensing potential 
(j) in harmonic space, treating these fields as Gaussian 
and computing power spectra using the publicly avail- 
able Boltzmann code CAMB [26 to ^^ax = 5000. Us- 
ing a fast spherical harmonic transform, we compute the 
unlensed temperature and polarization in pixel space, 
using an equicylindrical pixelization with {Ne^Nfj)) = 
(16384,32768) equally spaced points in (O^cj)). We then 
evaluate the lensed temperature and polarization fields 



X(n) =X(n + d(n)), 



(4) 



where X G {T, Q,^/} and tildes distinguish unlensed 
from lensed fields throughout, at each point of an N^ide = 
4096 Healpix pixelization [27], using cubic interpolation 
on the equicylindrical map to evaluate the right-hand 
side, and parallel translation of the spin-2 field Q ^ ill 
to transport polarization at the point (n + d(n)) to the 
point n. (We use an equicylindrical pixelization for the 



unlensed fields, rather than an irregular pixelization such 
as Healpix, so that interpolation is straightforward to im- 
plement.) Taking another spherical transform to obtain 
lensed T, and B maps in harmonic space, we compute 
lensed power spectra Cf^ for XY G {TT, TE, EE, BB}. 
As a memory optimization, we avoid storing full-sky 
maps by "striping" the sky into 16 latitude bands, and 
calculate the contribution to a^, a^, af^ from each 
band before moving onto the next. This allows each sim- 
ulation to fit onto a single core with ~2 GB memory. The 
above procedure is similar algorithmically to the publicly 
available code LensPix [28 , although the two codes differ 
in minor details of implementation. 

We first compare the mean over the N = 32768 real- 
izations 



XY 



1 

TV 



N 



EAXY 



(5) 



where XY e TT, TE, EE, BB to the predicted lensed 
power spectra computed by CAMB. With the resolution 
parameters given above, the lensed CMB power spec- 
tra of the simulations agree with CAMB's calculation of 
the lensed C^^'s to better than 0.1% for all spectra at 
Cax = 3000. 

We then estimate the covariance matrix between two 
different power spectra XY and ZW as 



1 ^ 

f^^^XY,ZW_ J- \^ AXY AZW nXY n 



ZW 



(6) 



Even in these noise- free simulations, the Gaussian ran- 
dom variance from the unlensed CMB makes the estimate 
of the covariance between individual multipoles noisy. 
We therefore further bin the power spectrum estimators 
into band powers 



XY 



XY 



(7) 



where Bf is a top hat function 



B 







i^<i< 

otherwise. 



(8) 



The band width is chosen to be sufficiently small so as 
to resolve the acoustic features in the spectrum. In prac- 
tice we take every multipole to ^ = 25 followed by uni- 
form bands of —li = 15. We choose not to bin the 25 
first multipoles as the derivatives of the power spectrum 
with respect to cosmological parameters exhibit strong 
variation at low multipoles, and averaging these varia- 
tions to one single bin at low-^ would give erroneous fi- 
nal results. The covariance matrix between these band 
estimators then becomes 



Gov; 



XY,ZW 



(9) 
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The Monte Carlo bandpower covariance is shown in 
Fig. [l] For visualization purposes, it is convenient to 
scale out the diagonal contributions by defining the cor- 
relation matrix 



R: 



XY,ZW 



Cov, 



XY.ZW 



Cov-. Cov^.^. 



(10) 



For display purposes, we also use a flat binning scheme 
by dividing the range of multipoles [2-3000] in 100 bins in 
Figs.[l}|5] As expected from previous studies, the covari- 
ance of the 5-modes is highly non-Gaussian [2T| [23l [29]. 
Interestingly, the EE^ BB and EE^ EE correlations in 
Fig. [l] are substantially larger than expected from the 
lowest-order analytic calculations in [HI |23] , and all but 
TT, TT show clear evidence for correlated structure on 
the acoustic scale that is again not expected. Although 
Ref. [21] also conducted simulation tests, their bands 
were much wider than the acoustic scale such that these 
structures were hidden. 



B. Analytic Approximation 

In order to develop a new analytic approximation to 
the covariance matrix, it is useful to first examine the 
BB^ BB correlation for which the existing models work 
well. The dominant terms in the analytic BB^ BB cor- 
relation expression can be compactly written as (cf. [29] 
Eq. 17) 



r^ BB,BB 




, (11) 



where unlensed CMB power spectra are denoted with 
tildes and Cf^ is the lensing potential power spectrum. 
Assuming that these fields are Gaussian, we can use the 
general prescription for Gaussian random fields G 



Gov 



GaGbyGcGd 



2e- 



GaGc ^GbG. 



^GaGd r^GbGc 



(12) 

for the unlensed CMB and (j) fields. 

To calculate power spectrum derivatives such as the 
ones appearing in Eq. (11), we take finite differences 



between lensed CMB power spectra computed using 
CAMB, rather than using a perturbative expansion in 
deflection angles. Since CAMB's algorithm for comput- 
ing lensed CMB power spectra includes terms of high or- 
der in deflection angles [30 , this approach to computing 
derivatives also includes high order terms, and in partic- 
ular does not break down at high I. Some implementa- 
tional details of the derivative calculation are presented 



in Appendix [Aj The model of Eq. ( 11 ) for the correlation 
matrix is shown in Fig. [2| 

Let us try to interpret the terms in Eq. ( [TT] ). The first 
term is the usual unconnected piece of the covariance 
that is the only term for a Gaussian random field. We 
will loosely refer to this term as the "Gaussian piece". 
The second and third terms involve the fact that the B 
field is constructed out of an unlensed E field and the 
lens potential field (j). In the second term, two BB band 
powers are connected by the covariance of the unlensed 
E fields they share. In the third term, they are connected 
by the shared (j) fields. Contributions to the correlation 
matrix for the second and third terms are shown sepa- 
rately in Fig. [3] 

The second term can therefore be interpreted as the 
covariance in BB band powers generated by cosmic vari- 
ance of the unlensed EE power spectrum. The covari- 
ance it generates is positive definite in that enhanced 
power in EE leads to enhanced BB across the spectrum 
thus correlating modes (see Fig. [sj right panel). 

The third term is the cosmic variance of the lens power. 
Here the correlation reflects the acoustic structure of the 
unlensed EE power spectrum. More power in the lenses 
allows more power from the acoustic peaks to transfer 
into 5-modes than the acoustic troughs (see Fig. [sj left 
panel). 



Finally, we note that Eq. (11) omits a fully connected 
term where the E and (j) fields are cross connected in- 
volving 4 unique multipoles rather than three. These 
contributions tend to sum incoherently and are subdom- 
inant in the covariance [23 . We omit this term in our 
analytic model. 

We can use these results to model the other co- 
variance terms. First consider BB^ XY where XY G 
TT, EE^ TE. In this case, there are no Gaussian or un- 
connected terms and 



Gov 



BB.XY 



dC, 



BB 

-=-^Cov 



EE,XY 



dC, 



XY 



dCf^ 



dcf 



-Gov? 



|. (13) 



In the perturbative limit for the deflection angles, this 
expression exactly models all terms in the covariance. 
However, again our expression has extended validity since 
the derivatives are evaluated nonlinear ly with CAMB. 

The case of XY = EE is illustrative as there is a sub- 
stantial correlation. The cosmic variance of the unlensed 
EE power spectrum produces contributions along the di- 
agonal but biased to a lower BB multipole ii < £2. This 
is due to the fact that most of the power in the low mul- 
tipoles of BB actually comes from where the unlensed 
EE spectrum peaks {ii ~ 1000). In previous analytic 
approaches, the term that was kept was for £2 = ^, which 
is linear in Cf^. 

Previous approaches have dropped the term associated 
with the cosmic variance of the lens power spectrum (the 
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FIG. 1. Monte Carlo covariance of the lensed CMB bandpowers computed from 32768 lensed CMB simulations in bands of 
Ai = 30. From left to right and top to bottom: TT, TE, EE, and BB. For visualization purposes we plot the correlation 
coefficient R defined in Eq. (10). The diagonal (of order unity) has been subtracted to enhance contrast. 



second term in Eq. (13)) under the justification that it 



is second order in cf^. In fact it is the dominant con- 
tribution to the covariance at ^1,^2 ^ 10^. This term 
causes a band structure in the EE dependence of the co- 
variance. Increasing the power in the lenses causes more 
power from acoustic peaks in EE to be transformed into 
BB power while also filling in power in EE at the acous- 
tic troughs. Thus peaks in EE are anticorrelated with 
BB and troughs are correlated. 

Finally, there are the cases for which XF, WZ G 
TT, TE^ EE. These cases are in principle more compli- 
cated in that even at the perturbative level, there are 
many terms that are not associated with the cosmic vari- 
ance of unlensed and lens potential power spectra. These 
are terms that connect the various unlensed, lensed and 
lens potential multipoles in the 4 point function. As in 
the case of BB^ BB we can again use the perturbative 



approximation as a guide. Here, there is a cancellation 
between the power spectrum covariance terms and the 
other terms associated with the unlensed fields for slowly 
varying unlensed power spectra. These other terms re- 
flect the fact that at high CMB multipole moment, the 
unlensed fields are all lensed by the same large scale lens 
realization. For a fixed lens, neighboring bands are anti- 
correlated by the exchange of power between them. This 
effect does not occur for the covariances with BB since 
there is no unlensed B field from which power can be 
taken. 



Given this close cancellation between terms associ- 
ated with the unlensed fields, we model only the cosmic 
variance of the lens power spectra in these cases. For 
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FIG. 2. Covariance of the lensed CMB bandpowers in bands of = 30, computed using the analytic model from §11 B| From 
left to right and top to bottom: TT, TE, EE, and BB. For visualization purposes we plot the correlation coefficient R defined 
in Eq. (10). Detailed comparison with the Monte Carlo covariance from Fig. [l] shows that the agreement is excellent. 



XY, WZ e TT, TE, EE 

XYWZ 1 



Cov^ 



21 



1 + J- 



dC, 



XY 



dCf 



-Govt 



.dC, 



WZ 



(14) 



In these cases the covariance takes a checkerboard pat- 
tern. For TT, TT or EE, EE enhanced lensing power 
makes modes near acoustic peaks smaller and larger near 
troughs. Thus peaks are correlated with peaks, troughs 
with troughs, and peaks are anticorrelated with troughs. 
Combining Eqs. (11), (13), (14), we have now devel- 



oped an analytic model for the lensed CMB bandpower 
covariance in all cases. Comparison with the Monte Carlo 



covariance from ^ 11 A shows that the difference is typi- 
cally < 10%, leading to discrepancies in parameter un- 



certainties on the order of 5% or less. We will explore 
this in more detail in §IV C[ 



III. MODEL INDEPENDENT LENSING 
INFORMATION 

As mentioned in the introduction, cosmological infor- 
mation from CMB lensing can be obtained either from 
precise measurements of lensed CMB power spectra, or 
by applying lens reconstruction techniques. In this sec- 
tion, we will quantify the relative amount of cosmolog- 
ical information which can be obtained using these two 
methods, in a model-independent way which uses Fisher 
information matrix techniques. 

In ^ III A[ we review the Fisher matrix formalism as 
applied to the lensing potential power spectrum. In ^IIIB 



Contributions of the individual terms in Eq. ( 11 ) to 

Left: cosmic 



FIG. 3 

the correlation matrix between BB bandpowers 
variance of the lens power (third term in Eq. (11)); right: 
cosmic variance of the unlensed EE spectrum (second term). 




FIG. 4. Rows of the correlation matrix defined in Eq. ( |1Q| ), 
between EE-EE band powers (top), and EE-BB (bottom), 
computed using either Monte Carlo simulations (solid lines) 
or our analytic model (dashed lines). Binning scheme follows 
Fig. ^ and the autocorrelation is omitted. 



we construct the Karhunen-Loeve (KL) basis to consider 
the relative information content. In §111 C[ we illustrate 
the necessity to take into account the non- Gaussian terms 
computed in ^TTj Finally in ^ III D , we apply our formalism 
to realistic CMB experiments. 



A. Fisher Information 

The Fisher information matrix quantifies the informa- 
tion in a given data set whose covariance matrix is known 
on a set of parameters Poc of interest. In order to quantify 
the lensing information in a model independent manner, 
instead of taking cosmological parameters we take the 
power spectrum Cf^ itself as the parameters of interest. 



1000 1500 



FIG. 5. Rows of the correlation matrix i?, as in Fig. [4] 
between BB-BB bandpowers computed using either Monte 
Carlo simulations (solid lines) or our analytic model (dashed 
lines). Binning scheme follows Fig. [l] and the autocorrelation 
is omitted. The upper curves are for i = 1306. 



The effect of any cosmological parameter of present or fu- 
ture interest can be thought of as a specific sum of these 
parameters. Rather than taking every ^ as a parameter, 
we follow [21 and implicitly assume that the power spec- 
trum is smooth in i so that we can approximate it with 
binned perturbations around the fiducial model. For each 
bin a in Cf^^ we define a parameter p^^ by 



In a 



£ Ifid ■ 



0<£< 



(15) 



where B^^^ describes the banding and is defined as 

" \ J otherwise. ^ ' 

In practice, we choose (j) bands with width Ai = 20 
(i.e. = 40 for £^^^ = 2000 and = 60 for 
^max = 3000). We define bands in InC/^ so that the 
power spectrum remains positive definite for large devi- 
ations. 

Note that any cosmological parameter variation that 
predicts a sufficiently smooth deviation from the fiducial 
model of S In Cf^ can be represented in these parameters 
as 

P« = 7T^y"<51nC/*B^'^ (17) 



(A£)« 



where (A^)^^ is the width of bin a. 

In general, given some data vector Dj which depends 
on parameters Pa , the Fisher matrix is given by 



oc(5 



ij 



dPi 

dpa 



(r^ D,D\- 
{CoVjj ) 



dDj 



dp ft ) ■ 



(18) 



We define a Fisher matrix by specializing to the case 
where the parameters p^ are the (j) bandpowers defined 
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in Eq. and the data vector Dj is the set of lensed 
CMB bandpowers Df^ (where XY G TT, TE, EE, BB) 
defined in Eq. ([7|). The derivatives {dDj /dpa) appear- 
ing in the Fisher matrix are computed nonhnearly using 
CAME, as described in the previous section. 

Unless otherwise stated, our model for the covari- 
ance matrix CovPj^ will be based on the semi-analytical 



model from [ II B (We show below that using the Monte 
Carlo covariance matrix from §11 A| gives essentially iden- 
tical parameter uncertainties.) We modify this all sky, 
cosmic variance limited covariance in two ways. First 
we include the possibility that the measurements contain 
Gaussian noise terms by replacing the Gaussian diagonal 
elements with 



Gov 



XY,X'Y' 



1 



2^ + 1 



where Nf-^ is the noise (cross) power spectrum. Second, 
we rescale the whole resulting matrix by l//sky where 
/sky is the fraction of sky covered by the data set under 
the usual assumption that the survey is sufficiently large 
that correlations across ~ 27r/^survey induced by the 
fundamental mode of the survey or sky cuts are irrele- 
vant across the acoustic separation. Also we are mainly 
interested in the effects of the non-Gaussian terms in the 
GMB covariance matrix at high I. An accurate forecast 
of parameter constraints or an analysis on data would re- 
quire a more detailed modeling of the shape of the survey 
{e.g. by introducing a minimum multipole ^min)- 

The Fisher matrix F^^ defined in this way quantifies 
the lensing information (in the form of constraints on 
the (j) bandpowers Pa) which can be obtained from noisy 
measurements of the lensed GMB power spectrum. 

We seek to compare this Fisher matrix to 



T^R 



E 



dCf 



dpa 



(Covg,)" 



dCf^ 



(20) 



the Fisher matrix of a direct reconstruction of . 

We will make the approximation that the covariance 
matrix Gov^/ of the reconstructed (j) bandpowers is given 
by the Gaussian expression 



Gov 



R 



/sky(2^+l) 



(C/^ + 7V/^)2(5,,,, (21) 



where Nf^ is the noise power spectrum of the reconstruc- 
tion. We use as a baseline two cases: a hypothetical cos- 
mic variance limited lens measurement where Nf^ = 0, 
and the idealized reconstruction noise coming from the 
lens reconstruction from quadratic combinations of GMB 
fields. For details on how this reconstruction noise is cal- 
culated, see 31 . 



The Gaussian approximation (21) makes several ap- 



proximations which we state explicitly. The power spec- 
trum of the lens reconstruction contains an off-diagonal 



contribution from (j) bandpowers [32] (this is the 'Wi" 
bias found by [33|) which should be folded into the co- 
variance Gov^/. It is also possible that there are contri- 
butions from higher-order terms in (j) [3T (this is the '^2" 
bias found by [35^); such contributions have been found to 
be small for temperature-based lens reconstruction, but 
this has not been checked for polarization. Finally, since 
the quadratic lens reconstruction is not a Gaussian field, 
its bandpower covariance may differ from the Gaussian 
expression (21). For temperature-based lens reconstruc- 
tion, this issue has been studied in [35] and the Gaussian 
expression has been found to be a good approximation 
(after slightly modifying the estimator along the lines of 
[36], see also [34] for alternate schemes), but the polar- 
ization case has not been studied. A complete treatment 
of these issues would be very interesting but is outside 
the scope of this paper; we will use the Gaussian approx- 
imation (21) as a first-order approximation to the exact 



Fisher matrix for lens reconstruction. 

Note that the inverse of the Fisher matrix is an ap- 
proximation for the covariance matrix of 



Gov, 



(F-i; 



al3 



(22) 



for both the band power (P) and reconstruction {R) 
Fisher matrices. 



B. Karhunen-Loeve Modes 

While the p^ basis of lens power spectrum perturba- 
tions is complete, it is not ideally suited for assessing the 
information content or analysis of data. Measurements 
of the many individual parameters would be highly noisy 
(and correlated, in the case where the lensed GMB band- 
powers Df^ are being used as the observable). In this 
section, we construct a more suitable basis whose eigen- 
modes are rank ordered in the relative information be- 
tween GMB bandpowers and lens reconstruction. More- 
over in the Fisher approximation, this basis provides a 
small set of relevant parameters whose errors are uncor- 
related for both bandpower and direct measurements. As 
a complete basis, it can be used to study any cosmological 
parameter which affects the lensing potential. 

To construct the eigenmodes, consider the Karhunen- 
Loeve (KL) transform: 



(23) 



where v^^^ and A*^*^ are the KL eigenvectors and eigenval- 
ues. We define the KL parameters rrik as linear combi- 
nations of the band perturbations or cosmological power 
spectrum deviations: 

mk = ^ v^^^Pa 

a 

= E-i'^(^E^i-^«" (24) 
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These KL modes have the property that their covari- 
ance, either measured from the lensed power spectra or 
from the reconstruction are diagonal and related by the 
KL eigenvalues. The eigenvectors are normalized such 
that all modes have unit variance for the direct recon- 
struction 



Gov 



Jkk' 



(25) 



and the KL construction then says that the eigenvalues 
are the relative variance from the CMB bandpowers 



C^^kk' 



kk' 



(26) 



The KL eigenvalues are therefore the ratio of the two 
covariances, and give a simple quantitative way to deter- 
mine which provides more information. If A^^^ < 1, the 
corresponding KL mode rrik is better constrained by the 
power spectra than by the direct reconstruction. 



TABLE I. Instrumental specifications used in this paper. Sen- 
sitivities are given in /xK arcmin. 



Name 


Frequency 


At 


Ap 


^FWHM 


/sky 


^max 




100 GHz 


81 


115 


9.5' 


0.8 


2000 


Planck 


143 GHz 


47 


79 


7.1' 


0.8 


2000 




217 GHz 


71 


122 


4.7' 


0.8 


2000 


Ground based 




1.0 


1.41 


1' 


0.016 


3000 


CMBPol 




1.0 


1.41 


1' 


0.8 


3000 



Amin = 1-11 when we use the covariance matrices com- 
puted from the simulations. We thus expect cosmological 
parameter errors to be modeled to better than a few per- 
cent for £max < 3000. We quantify this expectation for 



parameter examples in ^ IV C[ Our model captures the es- 
sential of the non-Gaussian structure of the lensed power 
spectra covariance. Most terms in this covariance model 
have been neglected in previous studies. 



C. CMB and Lens Cosmic Variance 



D. Finite Noise 



The KL construction allows a powerful test of physical 
self consistency of the bandpower covariance. If we con- 
sider an idealized measurement in which both the lensed 
GMB and the lensing potential are cosmic variance lim- 
ited (i.e. Nf^ = Nf^ = 0) then all KL eigenvalues must 
be > 1, since there cannot be more information in the 
lensed GMB bandpowers than the reconstruction. We 
also consider in this section that the full sky is observed, 

fsky — 1- 

If we treat the lensed B-mode as a Gaussian field 



(i.e. keep only the first term in the BB covariance (11)) 
then there are KL eigenvalues that strongly violate this 
physicality bound. For example, if we suppose that the 
BB power spectrum is measured to ^max = 2000, then 
we find Amin = 0.1. This problem disappears when we in- 
clude the full bandpower covariance: we find Amin = 1-4, 
showing that our covariance model passes this consis- 
tency test. These results are in agreement with p!2| |29]. 

At ^max = 3000, we find that all the polarization- 
related non-Gaussian covariances must be included in 
order to satisfy the physicality bound Amin ^ 1- For 
example, let us suppose that only EE and BB power 
spectra are measured (including TT and TE would only 
strengthen the example). If we make the Gaussian ap- 
proximation for Gov^^'^^, but use non-Gaussian values 
for Gov^^'^^ and Gov^^'^^, then we find Amin = 0.9 
and fail the consistency test. Analogously, if we make 
the Gaussian approximation for Gov^^'^^, but use non- 
Gaussian Cov^^'^^ and Gov^^'^^, then we find Amin = 
0.8 and fail. When we include the full non-Gaussian co- 
variance model from ^TTj then we do not find any viola- 
tion of physicality, even when all bandpowers TT, TE^ 
EE^ BB are included. In that case, considering all the 
covariance as Gaussian leads to Amin = 0.07. When 
the full non-Gaussian covariance from ^II] is used, we 
have Amin = 1.09 for the analytic covariance results, and 



While with a perfect reconstruction of the lensing po- 
tential power spectrum we cannot expect more informa- 
tion from the power spectra, considering a realistic re- 
construction with a finite noise could in principle lead to 
some modes which are better constrained by the power 
spectra than by reconstruction. 

In Table |lj we show some instrumental specifications 
that will be used throughout this paper. For Planck, we 
use the lowest three HFI frequencies with measured noise 
levels from [37 , with maximum multiple ^max = 2000 and 

/sky = 0.8. 

We also consider a futuristic CMB polarization satel- 
lite (denoted by "GMBpol") with a low noise level 
(At = 1/iK arcmin) and a resolution similar to the SPT- 
pol experiment [38 . Finally, we consider a ground exper- 
iment with instrumental characteristics from GMBpol as- 
suming a 650 deg^ survey and Planck sensitivity on the 
remainder of the Planck region (i.e. /sky = 0.784). We 
call this combination of Planck + ground-based experi- 
ment "Planck+G." 

The instrumental noise power spectrum for a single 
channel is [39j : 



N^^ = 



,/81n2 



(27) 



where XX = TT, EE, BB. For a multi-channel exper- 
iment, the noise power spectrum is Ni = 
where A^£(i) is the noise power spectrum of the i-th chan- 
nel. 

Note that even though GMBpol approaches the cosmic 
variance limit of the GMB, it does not reach the cosmic 
variance limit of lens reconstruction. The cosmic variance 
of the CMB fields themselves place an irreducible noise 
floor on even idealized reconstruction from quadratic es- 
timators. 
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FIG. 6. KL eigenvalues for Planck (solid/red) and CMBPol 



(dashed/green). Top: in 



2000, bottom: in 



: 3000. 



The KL eigenvalues are almost always larger than one, 
indicating that all the KL modes are better constrained 
by reconstruction than with the power spectra. The only 
exception is for CMBPol with a low cut-off at £max = 
2000. In that case Amin = 0.89, indicating that one KL 
mode is slightly better constrained by lensed CMB power 
spectra than by lens reconstruction. This number merely 
reflects the fact that given the low noise and beam of the 
CMBPol experiment, applying a cut-off at ^max = 2000 
degrades the ability of the lensing quadratic estimator to 
reconstruct the lensing potential. 

All the other KL eigenvalues are greater than one and 
they rapidly become much larger after the third mode 
(see Fig. [g]), indicating that only one or two KL modes 
actually contribute to the lensed CMB power spectra, 
which is in agreement with Ref. [21]. 



IV. CONSTRAINTS ON COSMOLOGICAL 
PARAMETERS 

Although the KL construction reveals extra informa- 
tion in lens reconstruction not available to lensed power 
spectra, accessing this information does not necessar- 
ily improve constraints on realistic cosmological parame- 
ters. Its impact depends on both how strongly and how 
uniquely cosmological parameter variations change the 
KL mode amplitudes corresponding to the new informa- 
tion. 

Our purpose is not to give exhaustive forecasts on cos- 
mological parameters for various experimental configu- 
rations. Rather, we wish to provide examples for when 
the extra KL information in reconstruction can and can- 
not make an impact. In § |IV A[ we define and test a 
means of comparing the two in the presence of parame- 
ters that change the acoustic peaks of the unlensed CMB. 
IV B[ we compare the errors from the power spec- 



In 

tra to the reconstruction and we assess the impact of the 
non- Gaussian covariance in ^ IIV C 



A. Additive Lensing Approach 

The KL mode decomposition is complete and hence 
the errors on the mode amplitudes rrik can be used to 
construct the Fisher matrices correspond- 
ing to lensed power spectra and lens reconstruction, for 
any set of cosmological parameters c: 

^ 1 drrik drrik 



KL,P 



^ A(^) dc dd ' 

k 



E 



drrik drrik 
dc dd 



(28) 



In other words, lensed CMB power spectrum constraints 
are downgraded (relative to lens reconstruction) by the 
KL eigenvalues. In practice, there are only a few eigen- 
values which are not ^ 1, and so lensed CMB power 
spectra are only sensitive to the first few eigenmodes. 

The derivatives {drrik /dc) appearing above can be 
computed from Eq. (|24|): 



drrik 
dc 



dliiCf 
~d^ 



(29) 



Some values of these derivatives with respect to neutrino 
mass, dark energy equation of state, and curvature are 
presented in Fig.j?! They show a general decreasing trend 
but are not monotonically decreasing. For example, the 
neutrino mass derivative with ^max = 2000 shows a sec- 
ond peak at the sixth KL mode (Fig. [7| left top panel). 

The KL Fisher matrices are designed to only account 
for the information carried by CMB lensing and should 
be added to any other source of information. Indeed 
these matrices are highly degenerate if the parameters 
that control the unlensed CMB acoustic peaks are al- 
lowed to vary. Two options are conceivable depending 
on the objective. If one wants to compare the ultimate 
amount of information on the lensing parameters carried 
by CMB lensing through the bandpower measurements 
or through a lens reconstruction, then fixing the high red- 
shift parameters in the KL matrices is a possibility. In 
this case, we expect the KL treatment to be fully accu- 
rate within the Fisher approximation, but the resultant 
error estimates are not meaningful unless other sources 
of information fix those parameters. 

The second approach is to add other sources of infor- 
mation to the Fisher matrix. The current leading source 
of information on these parameters is of course the acous- 
tic peaks themselves. The bulk of this information comes 
from the unlensed CMB spectra. To the extent that the 
lensing simply adds to the information in the unlensed 
CMB we can approximate the total Fisher matrix as the 
sum 

jpKL, PU ^ jpKL,P pU 



7U 



(30) 

where is the Fisher matrix constructed out of the un- 
lensed CMB fields with a Gaussian covariance. We call 
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FIG. 7. Derivatives of the KL modes ruk with respect to Q^h^ (top), w (middle), and Qk (bottom) with /max = 2000 (left 
column) and /max = 3000 (right column), for Planck (solid/red) and CMBpol (dashed/green). 



this the "additive lensing" approximation. For lensing 
parameters where lensing can actually destroy informa- 
tion in the power spectrum, such as Qk EI] and rrijy 
at sufficiently large values that the neutrinos are non- 
relativistic at recombination, this treatment is approxi- 
mate. On the other hand it is an approximation that 
affects the lensing reconstruction and power spectrum 
information alike. 

For the power spectrum information, there is a direct 
check of this approximation since we can construct the 
Fisher matrix from the lensed power spectra and the 
power spectrum covariance 



,dir,P 



IJ 



dPi 

dc 



(Gov 



D,D^-l 

IJ ) 



dPj 
dd 



(31) 



For the reconstruction, a direct check of the additive 
lensing approximation would require understanding the 
covariance between reconstruction and power spectrum 
statistics, as well as more subtle effects such as the A^i 
bias mentioned previously [40 . 

As an aside, note that for numerical stability when 
computing the power spectrum Fisher matrices it is im- 
portant to pick a parameter set c where the angular di- 
ameter distance degeneracy is manifest. Hence in prac- 
tice derivatives with respect to cosmological parameters 
are computed by adjusting the Hubble parameter h so 



that the acoustic scale is fixed when varying the values 
of other parameters. Our parameter basis is then com- 
posed of three lensing parameters i^rnj^^ and Qk) 
and six high-redshift parameters which control the un- 



lensed GMB: {Qch'^ .Qi.h'^ 



,0s}^ where Os is 



the angle subtended by the sound horizon at recombina- 
tion. These six high-redshift parameters are marginalized 
in all parameter constraints presented in this paper. 

We begin by testing the accuracy of the additive lens- 
ing approximation in the power spectrum case, where we 
can simply compare the Fisher matrix i^KL,Pt/ obtained 
in the additive lensing approximation (Eq. ([30|) to the 
exact Fisher matrix F^^^'-^ (Eq. (31)). As can be seen in 
Fig. |8j which presents the ratio of the errors computed by 
the two different techniques for a single additional lens 
parameter, the two approaches are not strictly equiva- 
lent. As expected, this is especially true for where 
the errors from the direct lensed Fisher matrix are typi- 
cally 5% larger than those predicted by our KL formalism 
and can approach 20% at high noise. For the differ- 
ence in the errors is constant at about 3% over the range 
of noise level considered. Finally, for ^ m^^ and the high 
fiducial value of 0.58eV, the agreement depends on the 
noise level. For very high levels of noise, most of the in- 
formation comes from the first few peaks of the unlensed 
GMB and that information can be reduced by lensing. 
For very low levels of noise, the lensing information sat- 
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urates to its sample variance level, while the unlensed 
CMB would have in principle retained information far 
out into the exponentially damped tail. In the interme- 
diate noise regime of interest to future CMB polarization 
experiments, the additive lensing approximation is ac- 
curate. Furthermore, we have explicitly verified that as 
the fiducial value for ^rriy is lowered, the discrepancy 
rapidly goes away. 

The same general trends apply to cases of multiple 
lensing parameters. As an example, we show in Fig. |9] 
(top panel), the constraints in the mj,) — w plane (VLk 
being fixed) for the Planck+G and CMBpol experiments. 
The ellipses from the two approaches are in very good 
agreement, thus validating the additive approach in the 
case where curvature is fixed. 

In the case where curvature is allowed to vary (bot- 
tom panel), parameters become highly degenerate in the 
lensed power spectrum effects and so the impact of un- 
lensed information becomes larger. For the CMBPol case 
the errors in (j{my) with w and curvature marginalized 
are larger by a factor of 1.26 and for the Planck+G exper- 
iment they are larger by 1.20 when comparing the exact 
to the additive approach. When making comparisons be- 
tween reconstruction and power spectrum information in 
such degenerate cases with curvature, one must bear in 
mind these curvature induced problems 



Eq. (30). This provides a metric for relative comparison 



of power spectrum and reconstruction lensing informa- 
tion. 




0.95 



0.90, 



4 6 

Aj, (/^K.arcmin) 



FIG. 8. Ratio of statistical errors computed using the addi- 
tive lensing approximation (i.e. ^^L,Pt7^ ^j^^ exact Fisher 
matrix (i.e. i^^^^^'^) for individual lensing parameters, as a 
function of noise level for a ^fwhm — l^^max = 3000 ex- 
periment. The error on each lensing parameter is computed 
with the other two lensing parameters fixed and high-redshift 
parameters marginalized. 

In summary, for the interesting cases where lensing 
provides most of the information on the lensing param- 
eters, the additive lensing approximation is accurate in 
the power spectrum case (i.e. the Fisher matrices i^K:L,PC/ 
and F^i^'^ agree). The additive lensing approximation is 
very convenient for comparing the cosmological informa- 
tion from lensed CMB power spectra and lens reconstruc- 
tion. Since the Fisher matrix is separated into a sum of 
unlensed and lensed contributions, we can simply com- 
pare the Fisher matrices i^KL,Pf/ pKL.RU ^Q^^ied in 




-0.9 



FIG. 9. 68% confidence limit (CL) ellipses in the {^ttIu) —w 
plane for the Planck+G (outer ellipses) and CMBpol (inner 
ellipses) experiments. Errors using the additive lensing ap- 



proximation (F 



are shown in dashed/black, and errors 



using the exact Fisher matrix (F'^^^'^) are shown in solid/red. 
Top: Qk is fixed. Bottom: Qk is marginalized. 



B. Power Spectra vs. Reconstruction 

Given the results in the last section, we can compare 
cosmological parameter constraints from lens reconstruc- 
tion and lensed CMB power spectra. Moreover, using the 
KL eigenmode formalism, we can explicitly verify how 
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FIG. 10. Statistical errors on ^rriiy (left), ^^(middle), and w (right) from lens reconstruction (solid/red) and lensed CMB 
power spectra (dashed/blue), using the first k KL eigenmodes (where I < k < 10), for CMBpol specs with ima^ = 3000. As 
expected, KL eigenmodes with /c > 3 do not contribute to the power spectrum constraints, but can contribute slightly in the 
lens reconstruction case. Statistical e rrors w ere computed using the "additive lensing" approximation and the Fisher matrices 

Top row: only one lensing parameter is varied and the other two are fixed; bottom 



F^L'^^ and F^L,i?t7 described in ^VA 



row: all three are varied with the two not shown marginalized. 



many KL modes actually carry the cosmological infor- 
mation. 

As a general statement, we find that reconstruction 
always carries more information than the lensed power 
spectra, regardless of the parameter considered. As the 
KL modes are rank-ordered in terms of highest relative 
information content in the power spectra, we can choose 
to truncate the summation defining the KL Fisher matri- 
ces (Eq. ([28])) to only use the information from the first 
few KL modes. Those cumulative errors are presented 
in Fig. [To] (upper) for cases where there is only a single 
additional lensing parameter. This case is the easiest to 
understand since the power spectrum information will be 
dominated by the first eigenmode. 

For CMBPol, the first KL eigenvalue is close to 1 and so 
the errors for any individual lensing parameter are com- 
parable. The power spectrum information saturates at 2- 
3 eigenmodes as expected. These are modes for which the 
extra information in the reconstruction becomes mani- 
fest, but for the three chosen lensing parameters the to- 
tal impact is small for and given the dominance 
of the first mode in the derivatives in Fig. ^ For the 
neutrinos there is a somewhat larger effect correspond- 
ing to large derivatives in both of the first two modes. 
Note also that for reconstruction, the eigenmodes are not 
rank-ordered so some of the higher modes can contribute 
more information than the lower modes. For neutrinos, 
the cumulative reduction of errors from the reconstruc- 
tion (relative to the power spectrum constraint) reaches 
0.7-0.8. 

When more than one lensing parameter is included, 
resulting degeneracies can make the information in re- 



construction more important. However the small deriva- 
tives in Fig. [T] still limits the practical relevance of this 
information, in that other sources like the unlensed CMB 
and more importantly other cosmological probes quickly 
dominate the net information. 

In Fig. [TT] we show a 2-dimensional example: com- 
parison of power spectrum and lens reconstruction con- 
straints in the (^ my)-w plane with VLk fixed, for Planck 
and CMBpol. Let us interpret this figure in light of our 
KL eigenmode construction. For Planck (upper panel), 
there is one (roughly vertical) direction which is con- 
strained by CMB lensing, and the lens reconstruction 
constraint is stronger than the lensed power spectrum 
constraint. There is another (roughly horizontal) direc- 
tion which is constrained by lens reconstruction, but very 
weakly constrained by lensing information in the power 
spectrum (it is constrained by the unlensed power spec- 
trum). This is consistent with the KL eigenvalues for 
Planck shown in Fig. [6j there is one KL eigenvalue which 
is a little larger than 1, and the second KL eigenvalue 
is 1. For CMBpol (bottom panel of Fig. 11), there is 
one direction where the lensed power spectrum and lens 
reconstruction constraints are nearly exactly equal, and 
another direction where lens reconstruction is somewhat 
better. This is consistent with the KL eigenvalues in 
Fig. [6j the lowest KL eigenvalue is almost exactly equal 
to 1, and the second KL eigenvalue is ~2. 

Similarly, although the KL analysis would imply that 
with 3 lensing parameters, there should be substantially 
better lensing reconstruction constraints, for the chosen 
parameters and their fiducial values the unlensed CMB 
information rapidly dominates. In Fig. 10 (lower), we 
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FIG. 11. Comparison of 68% CL ellipses in the (X^^zy)- 
w plane obtained from lens reconstruction (black/dashed) 
and lensed CMB power spectra (red/solid). Statistical er- 
rors are computed using the additive lensing approximation 
(i.e. Fisher matrices F^L,Ht7 p^i.,pu^ ^^^^ 

Planck (top panel) and CMBpol (bottom panel, ^max = 3000). 
The outer ellipse in the top panel is unlensed CMB constraint 
(not shown in the bottom panel since the unlensed constraint 
is much weaker than the lensed constraint). 



show the impact on the parameter errors of marginalizing 
the other two lensing parameters. In this case, all three 
types of lensing parameters show 10% or greater cumu- 
lative improvements from the reconstruction due to the 
higher modes breaking degeneracies, but they are still of 
the same order of magnitude as those of the power spec- 
tra. Further relative improvements here are limited by 
the unlensed CMB information which also weakly breaks 
these degeneracies in the additive approach. Note that 
we are somewhat underestimating the impact of the ex- 



tra reconstruction information when considering CMB- 
only sources of information, since this ability to break 
degeneracies in the unlensed CMB is degraded by lens- 
ing. Nonetheless the main point that in practice the ex- 
tra information accessible to lensing reconstruction with 
CMBPol is mainly orthogonal to realistic cosmological 
parameters remains. 



C. Impact of non-Gaussian covariance 

We conclude our analysis by investigating the impact 
of the non- Gaussian covariance of the lensed power spec- 
tra on the final errors on parameters. Since this does not 
involve reconstruction information, we work here with 
the direct lensed power spectra Fisher matrix i^dir,p^ 
rather than using the "additive lensing" approximation 
from 3lVA 



First note that if one considers the amplitude of the 
fiducial lensing spectrum Cf^ = A\QnsCf^\fid as the in- 
dependent lensing parameter as often done in the cur- 
rent literature [8^^, then non-Gaussian modeling is re- 
quired for any experiment that gains information from 
polarization. This follows from our KL treatment where 
the information on the amplitude comes almost exclu- 
sively from the first mode. For example for CMBPol and 
^max = 3000, (j{A\Qns) = 0.0011 for Gaussian covariance 
and a{A\Q^s) = 0.002395 for the non-Gaussian covariance 
from the simulations. Our analytic model captures this 
degradation to 2.8% yielding cr(Aiens) = 0.002329. 

If on the other hand one takes the lensing parameters 
as fundamental cosmological parameters, the impact of 
non-Gaussianity is hidden by marginalizing their impact 
on Cf^. In the previous studies where only the domi- 
nant BB-BB covariance was considered [21 the impact 
of non-Gaussianity on lensing parameters was small once 
the power spectra amplitude Ag and the dark matter 
density Vtch^ were marginalized for ^max < 2000. This 
simplification has been employed in the subsequent lit- 
erature to study parameter forecasts in a wider range of 
scenarios [42 . 

The analytic model allows us to separate out the im- 
pact of the new EE and TE covariance terms. In Fig.[T2j 
we show the net impact on the w — rrij, errors of in- 
cluding the non-Gaussian terms on the Planck+G (top) 
and CMBpol (bottom) experiments, with all the high 
redshift parameters marginalized but fixed. For the 
Planck+G experiment the overall impact is small as one 
would expect from just adding BB lensing information 
to Planck. For the CMBpol experiment, there is a more 
substantial effect. Interestingly this degradation is al- 
most entirely due to the new EE and TE terms and can 
be attributed to the use of lensing information in these 
spectra to break degeneracies between lensing parameters 
and high redshift parameters. 

These new terms can have an impact even if there is 
only one additional lensing parameter. For the CMBpol 
experiment if ^mjy and VLk are fixed, the error on w 
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CONCLUSIONS 
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FIG. 12. 68% CL ellipses in the {^mT^)-w plane, computed 
using the exact Fisher matrix F^^^'^ with fixed f^x, for the 
Planck+G (top) and CMBPol (bottom, ^^ax = 3000) experi- 
ments. Dot-dashed magenta lines are computed with the full 
Gaussian covariance matrix, solid red lines with the full non- 
Gaussian covariance from the semi-analytical model. Dashed 
black lines are computed using an intermediate covariance 
matrix which includes the non-Gaussian BB-BB covariance, 
but Gaussian covariance for all combinations of T and E. 



with the full non-Gaussian covariance is cr{w) = 0.118. 
However, if only the non-Gaussian BE — BB covariance 
is used, we have (t{w) = 0.095, a 24% difference. We 
conclude that if in the future lensing information out to 
^max = 3000 from the polarization fields becomes avail- 
able and dominates parameter errors, then all of the co- 
variance terms that involve polarization should be mod- 
eled for full accuracy. 



We have constructed a semi-analytic model of the co- 
variance matrix of the lensed power spectra of CMB tem- 
perature and polarization anisotropics. This model is 
able to reproduce the structure found in simulations of 
CMB lensing for the non- Gaussian terms in this covari- 
ance. More specifically, we have shown the existence and 
importance of second-order terms in the lensing poten- 
tial that were unaccounted for in previous studies [2T1 - 
|23] . Our model captures these effects and enables an effi- 
cient quantification of cosmological parameter errors that 
matches simulations to better than ~ 3% for ^max 

< 3000 

even in cases where the non-Gaussianity causes an order 
unity degradation in the errors. 

Using an parameter independent approach based on 
the decomposition of the information carried by CMB 
lensing in terms of Karhunen-Loeve eigenmodes, we have 
exhibited some cases where neglecting some of these sec- 
ond order terms lead to physical inconsistencies. These 
inconsistencies are removed once the covariance from our 
model is used. 

We then applied the KL eigenmode technique to com- 
pare the cosmological information that can be extracted 
either from measurements of the lensed power spectra or 
by reconstruction of the lensing potential using quadratic 
estimators. Although the non-Gaussian covariance of the 
lensed spectra has no significant impact on parameter 
errors for Planck, we found that it is non-negligible for 
forthcoming CMB experiments which will probe polar- 
ization at the arcminute scale. 

If the full non-Gaussian covariance is used then there is 
always more information, in principle, in the reconstruc- 
tion than in the lensed power spectra. In practice, the 
removal of various biases in the reconstruction and higher 
order terms in the reconstruction covariance matrix may 
degrade the final errors on parameters. Furthermore, this 
extra information is mainly in the detailed shape of the 
power spectrum of the lenses. Typical cosmological pa- 
rameters do not access this information as they mainly 
change the amplitude of the spectrum. 

The work presented here is one element in a joint and 
optimal likelihood analysis of CMB lensing. Our covari- 
ance model provides a computationally efficient means of 
calculating the covariance matrix of lensed CMB power 
spectra as a function of underlying cosmological or lens 
parameters. In the future, a full joint analysis will require 
more accurate techniques for a similar characterization 
of the reconstruction covariance as well as the covariance 
between power spectra and reconstruction observables. 
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On the other hand dC^^ /dCf^ cannot in general be 
approximated by a band response as it will have both 
a smooth piece from lensing and a d^^^i term from the 
unlensed CMB. However, to calculate the covariance it 
suffices to note that both the BB derivative and the Cov 
term are slowly varying on the ALq, scale. Thus the sum 
over £ e Ba means that we can replace the true derivative 
with the band average response again 



da 



XY 



da 



XY 



oa^ 



XY 



(A3) 



Appendix A: Derivative approximation for 
non- Gaussian covariance 



The covariance of BB with EE, TT or TE (see Eq.[T3 ) 
involves computing derivatives of lensed CMB power 
spectra with respect to unlensed spectra 



E 



BB 



da, 



EE 



-Gov 



f)r^y ' 

EE,XY ^^£2 



(Al) 



First note that daf_^^ / daf^ is a slowly varying function 
of both li and I given the broad kernel that transfers 
power between the EE and BB. Thus the derivative 
can be well approximated by the average response to an 
unlensed band perturbation of width ALq, = 10 for £ G 
band 



da 



BB 



With this replacement 



E 



dCf^ da 



XY 
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£eBa 



Cov 



EE.XY 



rEEr^XY 



,(A4) 



Note that while not necessary here, the derivatives can 
alternately be calculated more exactly by perturbing sin- 
gle ^'s on a sparse grid in the unlensed t. As a function 
of the unlensed t the derivatives can be separated as 



da, 



XY 



dar 



AfSf. 



Bp 



(A5) 



(A2) into two slowly varying pieces Ai and Bi 
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